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ABSTRACT 

We present accurate Monte-Carlo calculations of Lya radiation pressure in a range of 
models which represent galaxies during various epochs of our Universe. We show that 
the radiation force that Lya photons exert on hydrogen gas in the neutral intergalactic 
medium (IGM), that surrounds minihalos that host the first stars, may exceed gravity 
by orders of magnitude and drive supersonic winds. Lya radiation pressure may also 
dominate over gravity in the neutral IGM that surrounds the HII regions produced by 
the first galaxies. However, the radiation force is likely too weak to result in supersonic 
outflows in this case. Furthermore, we show that Lya radiation pressure may drive 
outflows in the interstellar medium of star forming galaxies that reach hundreds of 
km s~^. This mechanism could also operate at lower redshifts z <6, and may have 
already been indirectly detected in the spectral line shape of observed Lya emission 
lines. 
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1 INTRODUCTION 

HII regions around massive stars convert a significant frac- 
tion of the total bolo metric luminosity of young galaxies 
into Lya line emission ([Partridge fc Peebleslll967l : ISchaererl 
[2003). This Lya radiation can exert a large force on sur- 
rounding neutral gas, as the Lya transition has a cross- 
section that is ~ 7 orders of magnitude larger than the 
Thomson cross-section, when averaged over a freq uency 
band as wide as the resonance frequency itself (e.g. iLoebl 
I2OOII ). Not surprisingly, the impact of Lya radiation pres- 
sure on th e formation of galaxies has been di s cussed exten- 
sively (e.g[Cro3ll985l:lElitzur fc Ferlandlll986l: lBithelllll990l : 



lHaehnelt|]l995l : lOh fc HaimanI l2002l : iMcKee fc TanI 120081 ) . 
but the intricacies of Lya radiative transfer in 3D compli- 
cated an accurate numerical treatment of its dynamical ef- 
fect on the gas. Nevertheless, an approximate estimate can 
be obtained from simple energy considerations as shown be- 
low. 

Consider a self-gravitating gas cloud of total (baryons 
-|- dark matter) mass M and radius R that contains a cen- 
tral Lya source. The gravitational binding energy of the 
baryons inside the cloud, Eb ~ fii,GM^/(Om_R), can be 
compared to the total energy in the Lya radiation field in- 
side the cloud, Ea = La x ttrap- Here, La is the Lya lu- 
minosity of the central source (in erg s~^), and ttrap is the 
typical trapping time of Lya photons in the cloud owing to 
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scattering on hydrogen atoms. The Lya radiation pressure 
would unbind the baryonic gas fro m the cloud if Eg > Eb , 
i.e. La > Q.iGM'^KQ .n.Rttr:^^) (e.g. [Cffi3[l985l : lBithelllll9"90l : 
lOh fc HaimanI 12002! ) . In this approach, ftrap is one of the 
key paramet ers in se t ting t he Lya radiation pressure. Cal- 
culations by lAdamd (|l975l ) imply that itrap ~ 15fiight for 
3 ^ log TO ^5.5, and ttrap ^ 15(To/10^'^)^''^tiight otherwise, 
fo r a static, uniform, in finite slab of material (also see Fig 1 
of iBonilha et al.l 119791 ). Here to is the line center optical 
depth from the center to the edge of the slab, and iiight is 
the light crossing time if the medium were transparent (i.e. 
flight = R/c in the case of the cloud described above). Note 
however, that the precise value of itrap depends on other 
factors including for example, the gas distribution (dumpi- 
ness and geometry), the veloc ity distribution of th e gas, and 
the dust content of the cloud (jBonilha et al .11 19791 ). The Lya 
radiation pressure becomes comparable to gravity when 



= 1.1 






where La ~ La,ii x 10 erg s , and where we have 
substituted the virial radius of a galaxy mass M, 7?vir = 
0.97 kpc x(M/lO'*Af0)i/^(l + 2/16)~\ for 7? (Eq. 24 of 
Barkana & Loeb 2001). For comparison, a star forming 
galaxy can generate a Lya luminosity of La = (10''^ — 10*"^) x 
(SFR/M0 yr~^) erg s~^ , where the precise conversion factor 
depends on t he gas metallici ty and the stellar initial mass 
function (e.g. ISchaeren 120031 ). Therefore, a star formation 
rate of merely SFR ^0.01-0. IM0 yr~^ is needed to gener- 
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ate a Lya luminosity that is capable of unbinding gas from 
a halo of mass 10* - IO^'Mq. 

Halos of < 10^ Mq are very common at z ^6, and have 
a sufficiently large reservoir of baryons to sustain the above- 
mentioned star formation rates for a prolonged time. In this 
paper we provide a more detailed investigation of the mag- 
nitude of Lya radiation pressure in the environment of high- 
redshift star forming galaxies. In p articular, we use a L ya 
Monte-Carlo radiative transfer code (jDiikstra et al.ll2006r ) to 
compute Lya radiation pressure in a wider range of models. 
Our treatment of radiative transfer and our focus on the en- 
vironment of high-redshift star forming galaxies, distinguish 
this paper from previous work. We will show that the radia- 
tion force exerted by Lya photons on neutral hydrogen gas 
can exceed the gravitational force that binds the gas to its 
host galaxy by orders of magnitude, and may drive super- 
sonic outflows of neutral gas both in the intergalactic and 
the interstellar medium. 

The outline of this paper is as follows. In §[2] we describe 
how Lya radiation pressure is computed in the Monte-Carlo 
radiative transfer code, and show the tests that are per- 
formed to test the accuracy of the code. In § O we present 
our numerical results. Finally, § |4] summarizes the implica- 
tions of our work and our main conclusions. The cosmolog- 
ical parameter values used throughout our discussion are 
(Qm nA,nb,h) = (0.27,0.73,0.042,0.70) (JKomatsu et all 
l2008l ). 



where fcs is the Boltzmann constant, T the gas temperature, 
and rrip the proton mass. 

Under the assumption that /(jy, ^) has only a weak de- 
pendence on direction (which is reasonable given that Lya 
radiation scatters very frequently), I{u,yj) can be expressed 
as a first-order Taylor expansion in fx, i.e. I{v,fi) = a{y) + 
biy)^. In this so-called "Eddingt on approximation" , the ex- 
pression for flux simplifies to (e.g. lRvbicki fc Lightmanll979l . 
their Eq. 1.118) 

1 dj{v) c du[v 

3 dr ^ 



H{u) 



(4) 



127r dr ' 

where we have used the relation, u — AttJ/c, in which u{i') 
is the specific energy density in the radiation field at a fre- 
quency v. Furthermore, we have decoupled the gas' absorp- 
tion and emission functions from the Lya radiation field, 
and assumed that all neutral hydrogen atoms are in their 
electronic ground state (this assumption is justified in more 
detail in Appendix fX)) . i.e. n2 = and rii — uh- Under 
this assumption, dr = Xvdr — nH(T{v)dr with nu being the 
number density of neutral hydrogen atoms. Substituting this 
expression back into Eq. ((2} yields 



^ rsvd — 



1 d 
Suh dr 



dv u{v) 



1 dU 
3n_tf dr ' 



(5) 



where we defined U = J du u{v). Note that the cross-section 
does not appear in the final expression for the radiation 
forcfQ. 



2 LYq radiation PRESSURE 

The force Frad experienced by an atom in a direction n is 
related to the fiux through a plane normal to n. 



Fr, 



An 



du a[v)H{u), 



(2) 



where (j{v) is the Lya absorption cross-section at frequency 
V. The specific fiux is given by H{v) = ^ J dfi fil{fi,v), 
where I[v,ii) is the s pecific intensity of the radi ation field 
(see, e.g. Eq. 1.113 in lRvbicki fc Lightmanlll979r ). and /i = 
n • k in which k denotes the propagation direction of the ra- 
diation (i.e. /^ = 1 for radiation propagating perpendicular 
to the plane). 

The specific intensity obeys the radiative transfer equa- 
tion, which reads (in spherical coordinates) 



dl (1 - A*^) dl 



dr 



dfi 



X.{J ~ I) + S^r), 



(3) 



where in this equation fi = r ■ k/jr|, J{v) = J dfi I{^,i/) 
denotes the mean intensity, and Sv{r) the emission func- 
tion for newly created photons at frequency u and radius r 
(in photons cm^'^ s^^ sr^^ Hz"^, see e.g. Loeb & Rybicki 
1999). Furthermore, x^ = ^^77^ ./IaI , (^"■i - ^a) </'('") de- 
notes the opacity at frequency u (e.g. iRvbicki fc LightmanI 
1 19791 ). where ftp is Planck's constant, v^ = 2.46 x 10^^ Hz is 
the Lya frequency, ni(2) is the number density of hydrogen 
atoms in their electronic ground (first excited) state, B21 is 
the Einstein-B coefficient of the 2 — » 1 transition, (p{iy) is the 
line profile function (e.g. Rybicki fc Lightman 1979, their Eq. 
1.79), and Ai/a ~ ■^^. Here, Vth is the thermal velocity of 
the hydrogen atoms in the gas, given by vth = \/2k^3T/rrhp, 



2.1 Implementation in Monte-Carlo Technique 

In our Monte-Carlo simulation we sample the gas density 
and velocity fields with Ns = 5000 concentric spherical 
shells. The radius, thickness, and volume of shell j are de- 
noted by rj, drj, and Vj, respectively. We compute the ra- 
diation force using two approaches: 

• In the first approach, we calculate the energy density 
{U in Eq. O in the Lya radiation field as a function of 
radius: Using the Monte-Carlo simulation we compute the 
average time that photons spend in shell j, which we denote 
by {t)j. The total number of photons that is present in shell 
j at any given time is then given by Na,j — Na x {t)j, 
where A^o, is the rate at which photons are emitted. This 
yields the energy density, Uj = Najhua/Vj. Finally, Eq. ([SJ 
is used to compute the radiation force on atoms in shell 
j. Note that estimators of the energy density in -and the 
momentum transfer by- a radiation fiel d in a more gene ral 
context is discussed by e.g. ILucvI l|l999f ) and ILucvI ||2007| ). 

• In the second approach, we calculate the momentum 
transfer from a Lya photon to an atom in each scattering 
event, Ap = ftp(kin — kout)/27r. Here, kin and kout are the 
photons wavevectors before and after scattering. We com- 
pute the average total momentum transfer (i.e. summed over 
all scattering events) per photon in shell j, (AV)], and ob- 
tain the total momentum transfer from Pa.j = Na x {A'P)j. 
The force on an individual atom is obtained by divid- 
ing by the total number of hydrogen atoms in shell j, i.e 

^ The right-hand-side of Eq. JSJ is analogous to the usual pressure 
gradient force in fluid-dynamics which is not dependent on the 
scattering cross-section of the fluid particles. 
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Both methods should give identical results, provided 
that the Eddington approximation holds. 



2.2 Test Case: Sources in a Neutral Comoving 
IGM 

We begin by considering a Lya point source at a redshift 
z = 10 embedded in a neutral intergalactic medium (IGM) 
that is expanding with the Hubble flow. The photons scat- 
ter and diffuse away from the source while Hubble expansion 
redshifts the photons away from resonance. In this case, the 
angle- averaged intensity J iv) and its radial depe ndence can 
be calculated analytically l|Loeb fc Rvbickill999l ) . The avail- 
ability of analytic expressions for Jiy^r), and therefore the 
radiation force _Flad (through Eq.[5]), makes this a good test 
case for our code. 

In Figure [T] we plot the radial dependence of the energy 
density (in erg cm~'^) in the Lya radiation field for a model 
in which the central source is emitting Na,5i x 10^* photons 
s^^ (where we have introduced the dimensionless quantity 
iVa,54 = {Na/lOf'^ photons s~^). This corresponds to a lumi- 
nosity of La = A^a,54 X 1.6 X 10'*' ^ erg s~^, which re presents a 
bright Lya emitting galaxy fe.g. lOuchi et al.ll2008r ). The blue 
dotted line shows the energy density if the IGM were fully 
transparent to Lya radiation. In this hypothetical case all 
photons stream radially outward, and the energy density is 
given by La/(47rr^c). The red dashed line shows the energy 
density, U{r) — — J du J{v,r), derived from the analytic 
expression for ,J{r,v) given in Loeb & Rybicki (1999, their 
Eq. 21), while the black histogram shows the energy density 
extracted from the simulation (ij 12. 1[) . Clearly, the analytic 
and Monte-Carlo calculations yield consistent results. Scat- 
tering reduces the effective speed at which photons propa- 
gate radially outward, which enhances their energy density 
(especially at small radii) relative to the transparent case. At 
sufficiently large distances however, the photons have red- 
shifted far enough from resonance that they are propagat- 
ing almost freely to the observer, and the energy density ap- 
proaches La/^Tvr^c. We note that at a sufficiently high value 
of Na,54,, the fraction of hydrogen atoms that populate the 
2p (and 2s) levels is non-negligible and our assumption that 
(almost) all of the atoms populate their electronic ground 
state becomes invalid (so that the solution for U(r) in Fig- 
ure [T] breaks down). However, as we show in Appendix 1X1 
this only occurs when iVQ,,54 JjlO^, well beyond the regime 
considered in this paper. 

In Figure [5] we compare the radiation force to the grav- 
itational force on a single hydrogen atom, Fgrav = GM{< 
r)mp/r^, where M(< r) is the total mass enclosed within 
a radius r). We plot the ratio Frad/Jgrav scaled by M — 
IO^^MqQ The black dotted line (grey solid histogram) was 
calculated by applying Eq. ([5} to the energy density U{r) 
that was obtained by using the analytic (Monte-Carlo) ap- 
proach (also see Fig[T|. For comparison, the black solid his- 
togram was obtained by directly computing the momentum 
transfer rate from photons to atoms as outlined in § 12.11 



The number density of halos more massive than 10 Mq at 
2 = 10 is f^ lO"'^ comoving Mpc~^, implying that these rare 
halos are among the most massive ones in existence at that early 
cosmic time. 
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Figure 1. Radial profile of the energy density U{r) (erg cm^'') in 
the Lya radiation field surrounding a central source that is emit- 
ting Nci,iA X lO^** photons s~^ into an expanding neutral IGM. 
The blue dotted line shows C/(r) if the IGM were fully transpar- 
ent. The black solid (red dashed) line shows U{r) when radiative 
transfer is included using an analytic (Monte-Carlo) approach 
(see text). Scattering reduces the speed at which Lya photons 
are propagating radially outward, increasing U{r) relative to the 
transparent case. 
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Figure 2. The ratio of radiation to gravitational force on a hy- 
drogen atom as a function of physical radius in kpc. To scale 
out the dependence of this ratio on halo mass, M^, and produc- 
tion rate of Lya photons by the central source, Na, the vertical 
axis is normalized by Mf^ ii/Nc,,54 (see text). The black dotted 
(grey solid) line was obtained by applying Eq. JSJ to the energy 
density U{r) that was obtained by using the analytic (Monte- 
Carlo) approach (see Fig[TJ. The black solid line was obtained 
by directly computing the momentum transfer rate from photons 
to atoms in the Monte-Carlo code as outlined in § 12.11 The re- 
sults demonstrate that (i) the radiation force exceeds gravity at 
r < 10(7Va,54/M(j ii) kpc, and (ii) both methods yield consistent 
results. 



Figure [2] shows that the radiation force overwhelms gravity 
at small radii. The energy density scales approximately as 
91ogt//aiogr ~ -2.3 (FiglH). Therefore, F,ad/Fg,av oc r'^'^ 
and reaches unity at r ~ 10 physical kpc. 

The radiation force increases linearly with Nc, while 
the gravitational force scales as M. Thus, Frad/^grav scales 
linearly with the ratio 7?. = Na/M. To scale out the 
dependence on 7?., the vertical axis shows the quantity 
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(i^rad/Fgrav) X (Af n /iVc,54) , wheie Mn = (Af/1O"M0). For 
example, if Afn = 0.1 then Figure [2] shows that radiation 
pressure exceeds gravity out to r = 100 kpc, well beyond 
the virial radius of a halo of this mass at r^h ~ 6.6 kpc. 

Most importantly, Figure [2] shows that the two ap- 
proaches used to compute the radiation force in the simula- 
tion yield consistent results, with a noticeable deviation only 
at the largest radii (r ~ 1 Mpc). At large radii most photons 
stream outwards radially and the Eddington approximation 
that was used to derive Eq. (|4]) becomes increasingly unre- 
liable. 

Next, we use the radiative transfer code to explore the 
magnitude of the Lya radiation pressure for a range of mod- 
els which represent an evolutionary sequence of structure 
formation in the Universe. We focus on the Lya radiation 
pressure on gas surrounding (i) the first stars (^ I3.1|l : (ii) 
the first galaxies (§[3]2}; and {in) the interstellar medium of 
galaxies (§[23J. 



3 RESULTS 

3.1 Case I: A Single Massive Star in a Minihalo 

Numerical simulations of structure formation suggest the 
first stars that formed in our Universe were massive (Af* ~ 
lOOA^o), and formed as single objects in dark matter halos 
with masses of M ^ 10^ Mp, that cohapsed at z > 10 (e.i 
iHaiman et allll996l : lAbel et all I2OO2I : lYoshida et ai] I2OO6 
Here, we focus our attention on a star with a mass Af* = 
lOOAf© that formed at 2 = 15 in a dark matter mini-halo of 
mass M = 2 x lO'^Afg . The s tar emits 10^° ion izing photons 
per second (|Schaererl [20021 : I Abel et al.] 120071 ). We assume 
that the ionizing flux ionizes all the gas out to the virial 
radius of the dark matter halo (ryjr = .26 kpc) but not 
beyond that radius (IKitavama et al.ll2004l ). Hence, the IGM 
gas surrounding this central source is assumed to be neutral 
(iCHi = 1.0) and cold (Tgas = 300K, which corresponds to 
the temperature of the neutral IGM at « = 15 due to X-Ray 
heating, see e.g. Fig 1 of Pritchard & Loeb 2008). 

Recombination following photoionization converts ~ 
68% of all ionizing photons into Lya photons (|Osterbrockl 
[l989, p 387). Hence, the entire halo is a Lya source that is 
surrounded by neutral intergalactic gas. To determine the 
radial dependence of the Lya production rate {Si, in Eq.O, 
we need to specify the gas density profile. We assume that 
the gas distribution inside the dark matter halo is described 
by an NFW-profile with a concentration paramete r C = 5 
and a thermal corqj at r < 3rvir/4C (see Mailer fc Bullockl 
[2OOJ). We point out however, that our final results are not 
sensitive to our choice of S^(r). 

Once S„{r) has been determined, we find the radius, 
r, at which a Lya photon is generated in the Monte-Carlo 
simulation from the relation 



^-iV 



dr 47rr rinar, 



(6) 



^ With this gas density profile, the total recombination rate in- 
side the dark matter halo is Jq^" dr 47rr n^arcc ~ 4.5 X 10 
s~^ . The total recombination rate can be increased to balance 
the photoionization rate by introducing a clumping factor K = 
K>/(nH>2 ~ 2. 



where _R is a random number between and 1, TV = 
/„'^" dr 47rr^nHarcc is the total recombination rate inside 
the dark matter halo, and a^cc = 2.6 x 10~^^ cm'^ s~^ is the 
case-B recombination coeffi cient at a temperature T = 10'' 
K fe.g. iHui fc Gnedinlll997l ). Once the photon is generated, 
it scatters through the neutral IGM until it has redshifted 
far enough from resonance that it can escape to the observer. 

In the left panel of Figure [3] we show the energy den- 
sity (in erg cm~'^) of the Lya radiation field as a function 
of radius. The red solid line represents a model in which we 
assumed the IGM to follow the mean density and Hubble ex- 
pansion right outside the virial radius. The blue dotted line 
shows a more realistic model in which the IGM is still over- 
dense near the virial radius, and in which the intergalactic 
gas is gravitationally pulled towards the minihalo (see Dijk- 
stra et al 2007 for a quantitative description of the density 
and velocity profiles based on the model of Barkana 2004). 
The black dashed line shows the same model as the red solid 
Ime but with the neutral fraction increasing linearly between 
rvir and 2rvir. This provides a better representation of the 
fact that the central population HI star emits ionizing pho- 
tons with energies ^54 eV, which can photoionize hydrogen 
(and helium) atoms that lie deeper in the IGM. The goal of 
this model is to investigate whether our results depend sen- 
sitively on the presence of a sharp boundary between HI and 
HII. 

All models show that the radiation energy density 
within the fully ionized minihalo (r ;^rvir = 0.26 kpc) has 
only a weak dependence on radius, i.e. dlog U/dlogr ^ — 1. 
Naively, this may appear surprising given the fact that 
within the model, no scattering occurs within the virial ra- 
dius and one may expect the energy density in the Lya 
radiation field to scale as f/ oc r~^. However, in reality U 
obtains only a weak radial dependence because the radia- 
tion can be scattered back into the ionized minihalo as soon 
as it 'hits' the wall of neutral IGM gas. Lya photons are 
therefore trapped inside the ionized minihalo and their en- 
ergy density is boosted to a value that is only weakly de- 
pendent on radius. On the other hand, for r ^^vir we find 
that dlogU/dlogr <. — 2, which is because Lya photons 
are trapped more efficiently near the edge of the HII re- 
gion, while they stream freely outwards at larger radii (as 
in § 12.21 and Fig [l| . Figure [3] shows clearly that the radial 
dependence of the Lya energy density is not sensitive to the 
detailed model assumptions about the gas in the IGM. 

In the right panel of Figure |3] we show the ratio be- 
tween the radiation force (Eq. [5} and the gravitational force 



on a single hydrogen atom) , _Fg, 



GM{< r)mp/r , where 



A/(< r) is the total (baryons + dark matter) mass enclosed 
within a radius r. In all models, radiation pressure domi- 
nates over gravity by as much as ^2 orders of magnitude. 
The radiation force is largest for the models in which the 
IGM is assumed to be at mean density, because of the nj^ 
factor in the equation for the radiation force (Eq. [S|. Note 
that the spike near r ~ 2.6 kpc for the other two models is 
due to an artificial discontinuity in the IGM velocity field 
that exists in this model. 

Lya radiation pressure may operate throughout the life- 
time of th e central star. Over a lifetime of ~ 2.5 Myr (see 
Table 4 of lSchaererll2002l ). this mechanism is capable of ac- 
celerating the gas to velocities of 10 (50) km s~^ at r = 0.3 
kpc in the model represented by the blue dotted {red solid) 
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Figure 3. The energy density in tlie Lya radiation field (left panel), and tlie ratio between tlie radiation and the gravitational forces 
{right panel) for a single very massive star (M, = IOOMq) in a minihalo Mj, = 2 X 10® Mq . The red solid line represents a model 
in which the IGM is at the mean cosmic density and undergoes Hubble expansion right outside the virial radius at r = 0.26 kpc. The 
blue dotted line shows a more realistic model in which the IGM is overdense near the virial radius, and in which the intergalactic gas 
is gravitationally pulled towards the minihalo (see text). The black dashed line shows the same model as the red solid line but with 
the neutral fraction increasing linearly between r^iT and 2rvir. The star ionizes all the gas out to the virial radius. Lya photons freely 
propagate until they reach the edge of the HII region, where they are likely to be scattered back into the ionized minihalo. This yields 
an almost constant radiation energy density. Once outside the HII region, the radiation force dominates over gravity out to r = 10 kpc 
and may accelerate neutral gas outside the HII region to velocities of order ~ 10 km s~^. 
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Figure 4. Same as Figure [S] but for the case of a star forming galaxy (M* = 0.34Mq yr~^, see text), surrounded by an HII region 
with a radius -Rhii = 50 kpc (Rhii = 20 kpc) for the solid line (dotted line), which is in turn surrounded by a fully neutral intergalactic 
medium (IGM). For the assumed total halo mass of Mtot = 10^ Mq, the pressure exerted by the Lya photons is not large enough to 
exceed gravity. However, radiation pressure wins for Mtot = 10* M©, but even in this case the radiation force is not large enough to 
produce a significant wind speed in the IGM (see text). 



line, and to 4 (16) km s~^ at r = 0.4 kpc (the reason for 
this large difference is that the edge of the HII region hes at 
r — 0.26 kpc. Hence, gas at r = 0.3 kpc is separated by 0.04 
kpc from this edge, while gas at r = 0.4 kpc is separated by 
a distance that is 3.5 times larger). 

Thus, Lya radiation pressure can accelerate the gas to 
velocities that exceed the escape velocity from the dark mat- 
ter halo (wesc ~ \/2ucirc ~ 8 km s~^) as well as the sound 
speed of the intergalactic medium (c^ = 2.2(T'gas/300 K)^^^ 
kms-i). 

Note the as the gas is pushed out and its velocity pro- 
file changes, the subsequent radiative transfer is altered. For 
example, we repeated the radiative transfer calculation for 
models in which gas at r^ir < r <;2rvir was accelerated to 
velocities in the range 10 — 20 km s^^ (outward) and found 
a slightly shallower profile for U{r) which lowered the radi- 



ation force by a factor of ^ 3. Consequently, the accelera- 
tion of the gas decreases with time and the actual velocities 
reached by the gas are lower than the estimates given above 
by a factor of a few. Nevertheless, the resulting velocities 
are still substantial. 

Our calculations imply that Lya radiation pressure can 
affect the gas dynamics in the IGM surrounding minihalos 
that contain the first stars. The impact of Lya radiation 
pressure increases with decreasing density of the surround- 
ing gas in the IGM. In practice, the distribution of the IGM 
is not spherically symmetric. Instead, the density is expected 
to vary from sightline to sightline (being large along fila- 
ments and small along voids). Our results imply that Lya 
radiation pressure will be most efficient in 'blowing out' the 
lower density gas. This conjecture is supported by the ten- 
dency of Lya photons to preferentially scatter through the 
low-density gas; their propagation along the path of least re- 
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sistance would naturally boost up the Lya flux there. This 
effect will be moderated by the tendency of the HIT region 
around the first stars to extend further into the low density 
gas (in 'butterfly'-like patterns, e.g. Abel et al, 1999). 

If the central star dies in a supernova explosion, then 
the resulting violent outflow could blow most of the baryons 
out from the minihalo. However, stars with masses in the 
range SOA'/© <.A4, ^14OM0 and M* ^26OM0, are not ex- 
pected to end their lives in a super nova. Instead, these stars 
collapse directly to a black hole (JHeeer fc WooslevI |2003) 
and have weak winds (because of the lack of heavy elements 
in their atmosphere), so that radiation pressure may be the 
dominant process that affects their surrounding IGM. 

In summary, Lya radiation pressure on the neutral IGM 
around minihalos in which the first stars form, can exceed 
gravity by orders of magnitude and launch supersonic winds. 
Our limited analysis does not allow a detailed discussion 
on the consequences of these winds. This requires 3D sim- 
ulations with cosmological initial conditions that capture 
the full IGM density field around the minihalo and that 
track the evolution of the shocks that may form in the IGM. 
Such simulation are numerically challenging as they require 
self-consistent treatment of gas dynamics and Lya radiative 
transfer in a moving inhomogeneous medium. 



3.2 Case II: A Young Star Forming Galaxy 

Our second case concerns a young galaxy that is forming 
multiple stars in a dark matter halo of mass M = IO^Mq 
at z = 10. We assume that the galaxy is converting a 
fraction f, = 10% of its baryons into stars over ~ O.ltu 
l|Wvithe fc Loebll2006l ). where tn = [2/3H{z)] ~ 0.49 Gyr, 
is the age of the Universe at z = 10. This translates to a 
star formation rate of M* = 0.34Mq yr~^. For population 
III stars forming out of pristine gas, the t otal emission ra te 
of ionizing photons is iVion ~ 3 x 10^^ s"^ (|Schaereill200^ fl. 
If ~ 1% of the ionizing photon s escape from the galaxy 
IjChen et al.ll2007l : iGnedin et ail i2008'l. then this translates 
to a Lya luminosity of La = 3 x 10*^ erg s"'^. Furthermore, 
this galaxy can photoionize a spherical HII region of a radius 
^Hii ~ 50 physical kpc. Note however, that other ionizing 
sources would likely exist within this HII region. Indeed, 
clusters of sources are thought to determine the growth of 
ionized bubbles during reionization. This results in a charac- 
teristic HII region size that is significantly larger than that 
produced by single s ource, especially during the later stages 
of re ionization (e.g. iFurlanetto et al.l |2004| : iMcQuinn et al.l 
120071 ). In this framework, our model represents a star form- 
ing galaxy during the early stages of reionization or alter- 
natively a galaxy that lies 50 kpc away from the edge of a 
larger ionized bubble. 

In this particular case, the majority of all recombination 
events occur in the central galaxy. Thus, we initiate all Lya 
photons at r = in the Monte-Carlo simulation. We assume 
that the gas is completely ionized out to 7?Hn ~ 50 kpc, 
beyond which it is neutral. As shown in § 13.11 this abrupt 



* More precisely, the ionizing photon production rate is ~ 10^* 
s~^ x(SFR/M0 yr~^) in the no-mass-loss model of Schaerer 
(2002) in which metal- free stars form according to a Salpeter IMF 
with Miow = IMq and Afhigh = 500M@ (his model 'B'). 



transition in the ionized fraction of H in the gas does not 
affect our results. 

The left panel of Figure |4] shows the energy density (in 
erg cm~'^) of the Lya radiation field as a function of radius. 
The solid line represents the model discussed above. A kink 
in the energy density is seen at the edge of the HII region 
(see §[3TT]for a more detailed discussion of the profile). The 
dotted line represents a variant of the model in which we 
have reduced the size of the HII region to -Rhii ~ 20 kpc. 

The right panel of Figure |4] shows the ratio between the 
radiation and the gravitational forces on a single hydrogen 
atom. In our fiducial model, the radiation force does not 
exceed gravity; rather, at the edge of the HII region, gravity 
is ~ 3 times stronger. The radiation force becomes equal to 
the gravitational force if -Rhii ~ 20 kpc. This requires an 
extremely low [by a factor ~ (50/20)^] escape fraction of 
ionizing photons, /osc ~ 6 x 10~*. 

Alternatively, radiation pressure is important when the 
halo mass of the star forming region is reduced to 10^ Mq. 
Halos of this mass are the the most abundant halos at z ~ 10 
that are capable of cooling via excitation of atomic hydro- 
gen (i.e. their virial temperature just exceeds Tvir ~ 10'* K, 
e.g. Barkana & Loeb 2001). The total gas reservoir inside 
these halos is Afb = n^-^tot ~ 1.5 x IO^'Mq, and so these 
halos can sustain a star formation rate of M, — O.3A/0 yr~^ 
for up to ~ 50 Myr. However, even if the radiation force is 
allowed to operate for ~ 50 Myr, we find that radiation pres- 
sure cannot accelerate the gas in the IGM to velocities that 
exceed ~ 1 km s~*. We therefore conclude that although 
Lya pressure may exceed gravity in the neutral IGM that 
surrounds HII regions around Mtot ~ 10^ Mq halos, the ab- 
solute magnitude of the radiation force is too weak to drive 
the IGM to supersonic velocities. 



3.3 Case III: Lya Driven Galactic Supershells 

In principle, Lya radiation pressure can be important when 
neutral gas exists in close proximity to a luminous Lya 
source. So far, we focused our attention on HI gas in the 
IGM. However, neutral gas in the interstellar medium (ISM) 
of the host galaxy is located closer to the Lya sources and 
should be exposed to an even stronger Lya rad iation pres- 
sure. Indeed, it has been d emonstrated (e.g. lAhn fc Led 
|2002| : IVerhamme et al.l l200a ) that scattering of Lya pho- 
tons by neutral hydrogen atoms in a thin (with a thickness 
much smaller than its radius), outflowing 'supershell' of HI 
gas surrounding the star forming regions can naturally ex- 
plain two observed phenomena: (i) the common shift of the 
Lya emission line towards the red relative to metal absorp- 
tion lines and the host galaxy's systemic redshi ft determined 
from other nebular recombination lines (e.g. IPettini et al.l 
2001; Shaplcy ct al. 2003); and (li) the asymmetry of the 
Lya line with emission extending well in to its red wing (e.g. 
iLeaueuxetalll 19951 : iTapken et aLlbOOTh . 

The existence of thin, outflowing shells of neutral 
atomic hydrogen around HII regions is confirmed by HI- 
observations of our ow n Milky- Way (Heile s 1984) and other 
nearby galaxies (e.g. iRvder et al.l Il995l ). The largest of 
these she Us, so-called 'supe r shells', have radii of rmax ~ 1 
kpc (e.g iRvder et all 1 19951 : iMcClure-Griffiths et all I2002D 
and HI column densities in the range TVhi ~ 10*^- 
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Figure 5. Same as in Figure |3] and Figure |4l but for models in which a central Lya source of luminosity L^ is surrounded by a thin 
(''sh min = 0-9''sh max) Spherical shell of HI that is expanding at Dgj^ = 200 km s^^. The blue dotted lines {red solid lines) represent 
a model in which A^hi = 10^^ cm~^ and rgj, ^^x = 0.1 kpc (A^hi = 10^® cm~^ and rsj, m^x = 1-0 kpc). When calculating Fgrav, we 
assumed the total mass enclosed by the supershell to be 10* M© (see text). The results strongly suggest that Lyo radiation pressure 
may be dynamically important in the interstellar medium of galaxies. The total radiation force on the shell (obtained as a sum over all 
atoms) may be computed via _Ftot = MpLa/c, where Mp may be thought of as a force multiplication factor that depends both on the 
shell's outflow speed, rioxp, and its HI column density, A'^hi- The dependence of Mp on these parameters is shown in Fig|6] 



10^^ cm"^ (e.g iLequeux et al.1 1 19951 : iKunth et all 1 19981 : 
IVerhamme et al.ll2008i ). Supershells are thought to be gen- 
erated by stellar winds or supernovae explosions which 
swe ep-up gas into a thin expanding neutral shell (see 
e.g. IXenorio-Tagle fc Bo denheimer 19881, for a review). The 
back-scattering mechanism attributes both the redshift and 
asymmetry of the Lya line to the Doppler boost that 
Lya photons undergo as they scatter off the outflow on 
the far side of the galaxy back towards the observer (e.g. 
Lee fc AhnI Il998l: lAhii fc Leel |2002| : lAhn et al.l l2003l : lAhnl 
20041 : IVerhamme et al.ll2006l . mom . It is interesting to in- 



vestigate whether Lya radiation pressure may provide an 
alternative mechanism that determines the supershell kine- 
matics. 

In Figure [5] we show the energy density {left panel) and 
the Lya radiation force {right panel) for two models. Both 
models assume that: (i) there is a Lya source at r = 
with a luminosity of La = 10*'' erg s~^; {ii) the emitted 
Lya spectrum prior to scattering has a Gaussian shape as a 
function of photon frequency with a Doppler velocity width 
of (T = 50 km s^^; {Hi) the spatial width of the supershell is 
10% of its radius; and (iv) the shell has an outflow velocity of 
V = 200 km s~^. The blue dotted (red solid) lines represent 
a model in which the supershell has a column density of 
A^Hi ~ 10^^ (A'^Hi ~ 10^®) cm~^ and a maximum radius that 
is Tsh = 0.1 kpc (rah = 1.0) kpc. Our calculations assume 
that there is no neutral gas (or dust) interior to the HI 
supershell. 

The left panel of Figure [5] shows that inside the super- 
shell the energy density decreases more gradually than r~^ 
because of photon trapping (similarly to the previously dis- 
cussed cases in §[3TT]-§[3]2]). The shell with the larger column 
of HI is more efficient at trapping the Lya photons, and thus 
yields a flatter energy density profile. In both models the en- 
ergy density drops steeply within the supershell (the energy 
density decreases as r~^ outside the shell, if no scattering 
occurs here). 

The right panel of Figure [5] shows the ratio between 
the radiation and gravitational forces. Towards the center 
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Figure 6. The multiplication factor Mp provides the total force 
that Lya photons exert on a spherical HI shell, _Ftot = MpLa/c, 
where La is the Lya luminosity of the source in erg s~^ and c is 
the speed of light. The plot shows Mp as a function of the expan- 
sion velocity of the HI shell, Ugij, for three values of HI column 
density, and under the assumption that the HI shell surrounds 
an empty cavity. The parameter Mp provides a measure of the 
efficiency by which Lyo photons can be 'trapped' by the shell of 
HI gas. Thus, Mp increases with increasing Afjji and decreasing 
v^l^. The dotted horizontal lines show the values of Mp that have 
been derived in the past (for a static, uniform, infinite slab of 
material e.g. Adams 1975). 



of the dark matter halo, baryons dominate the mass density 
and an evaluation of Fgrav requires assumptions about the 
radial distribution of the baryons. For simplicity, we consider 
a fixed total mass interior to the supershell of M(< r) = 
10** Mq, so that Fgrav(r) = GM{< r)mp/r^. Note that any 
assumed mass profile M{< r) will not affect the results as 

long as Frad > -Fgrav . 

We find that the radiation force exceeds gravity in 



both examples under consideration. For A'^hi 



lo-" 
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and Tsh ~ 0.1 kpc, F,.^ 



10-Fgrav Thus, radiation pres- 



sure would have been important even if we had chosen 
M ~ 10^ Mq. In the model with Nm 



rah = 1.0 kpc, Fy, 



10-Fgrav, and radiation pressure would 



have been important even if M ~ 10 Mq. Hence, our cal- 
culations strongly suggest that Lya radiation pressure may 
be dynamically important in the ISM of galaxies. 

The total Lya radiation force is obtained by summing 
the force over all atoms in the supershell. It is interesting 
to compare this force to La/c. The latter quantity denotes 
the total momentum transfer rate (force) from the Lya ra- 
diation field to the supershell under the assumption that 
each Lya photon is re-emitted isotropically after entering 
the shell (including multiple scatterings inside the super- 
shell). In Figure [6] we plot the quantity Mf which is defined 



Mf 



Ea 



La/c 



(7) 



as a function of the expansion velocity of the shell, v^h for 
three different values of A'^hi. 

Figure |6] shows that Mf , which can be thought of as 
a force multiplication factoqj, greatly exceeds unity for low 
shell velocities and large HI column densities. The parameter 
AIf is related to the mean number of times that a Lya pho- 
ton 'bounces' back and forth between opposite sides of the 
expanding shell. For example, Mf ~ 1 when all Lya photons 
enter the shell, scatter once, and then escape from the shell 
in no preferred direction. On the other hand, Mf = 3 when 
all Lya photons enter the shell, scatter back towards the op- 
posite direction, and then escape in no preferred direction 
after scattering in the shell for a second time. A schematic 
illustration of this argument is provided in Figure [7] Note 
that when Mf = 1 {Mf ~ 3), each photon spends on aver- 
age a timescale of rsh/c (3rsh/c) in the bubble enclosed by 
the shell. In other words, the factor Mf relates to the 'trap- 
ping time', ttrap, that denotes the total time over which Lya 
photons are trapped inside the supershelO (see § [1} through 
the relation Mf = itrap/(rsh/c). Indeed, when Vah — > we 
find that Mf reproduces the value 15(ro/10 ^'^) (indicated 
by horizontal dotted lines ) that was found bv lAdamg (| 19751 ') 
and iBonilha et al.l (|l979l ) reasonably well (keeping in mind 



^ This term derives from the (time-depen dent) force- 
multip lication function M{t) that was introduced bv lCastor et al.l 
1I1975I V as -Frad = Af (t)(reLboi/c). Here, F^ad is the total force 
that radiation exerts on a medium, Te is the total optical 
depth to electron scattering through this medium. The function 
M{t) arises because of the contribution of numerous metal 
absorption lines to the medium's opacity, and ca n be as large 
as M max(f) ^ 10'^ in the atmospheres of O-stars 11 Castor et al.l 
I1975II . 

^ This argument ignores the time spent on scattering inside the 
supershell itself. Photons penetrate on average an optical depth 
T = 1 into the shell. If this corresponds to a physical distance that 
is significantly smaller than the thickness of the shell (denoted by 
Argil ) ! then only a tiny fraction of the photons will diffuse through 
the shell. Hence, when averaged over these photons, ignoring the 
time spent inside the supershell itself is justified. Alternatively, 
photons with a mean free path that is at least comparable to 
the thickness of the shell, only spend a time ~ Argh/c <C r^h/c 
inside the supershell, which provides a negligible contribution to 
the trapping time. 




Mp^=l 



Mf=3 



Figure 7. A schematic illustration of the origin of the force mul- 
tiplication factor Mf in an expanding supershell. The Lya source 
is located at the center of the expanding HI supershell. The solid 
line represents the trajectories of photons that enter the shell, 
scatter once, and then escape from the shell isotropically (indi- 
cated by the dashed lines). This corresponds to Mf = 1. On 
the other hand, the dot-dashed line represents the trajectories of 
photons that enter the shell, are scattered back in the opposite 
direction, and escape with no preferred direction after scattering 
in the shell for a second time. This corresponds to a case with 
Mf = 3. In general, Mf = ttrap/(rsh/c). 



that these authors derived their result for a static, uniform, 
infinite slab of material, and assumed different frequency 
distributions for the emitted Lya photons). 



4 CONCLUSIONS 

We have applied an existing Monte-Carlo Lya radia- 
tive transfer code (described and tested extensively in 
iDiikstra et al.ll2006i ) to the calculation of the pressure that is 
exerted by Lya photons on an optically thick medium. This 
code enabled us to perform (the first) direct, accurate cal- 
culations of Lya radiation pressure, which distinguishes this 
work from previous discussions on the importance of Lya ra- 
diation pressure in various astrophysical environments. 

We have focused on a range of models which repre- 
sent galaxies at different cosmological epochs. In § 13.11 we 
have shown that the Lya radiation pressure exerted on the 
neutral intergalactic medium (IGM) surrounding minihalos 
(Aftot ~ 10'' Af©) in which the first stars form, can exceed 
gravity by 2-3 orders of magnitude (Fig [3)l , and in princi- 
ple accelerate the gas in the IGM to tens of km s"'^. Thus, 
Lya radiation pressure can launch supersonic winds in the 
IGM surrounding the first stars. Our analysis did not allow 
a detailed study of the consequences of these winds. A com- 
prehensive study would require numerical simulations that 
capture the full IGM density field around minihalos in 3D 
and track the evolution of the shocks that may form in the 
IGM together with the Lya radiative transfer. In this paper, 
we have also shown that Lya radiation pressure is important 
in the neutral IGM that surrounds the HII regions produced 
by galaxies with a total halo mass of Aftot = 10* Af© (FigU) 
These are the lowest mass, and hence the most abundant, 
halos in which gas can cool via atomic line excitation. Here, 
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however, the absolute magnitude of the radiation force is 
too weak to drive the gas to supersonic velocities. 

Finally, we have shown in § 13.31 that the Lya radiation 
pressure exerted on neutral gas in the interstellar medium 
(ISM) of a galaxy can also have strong dynamical conse- 
quences. In particular, we have found that the Lya radiation 
force exerted on an expanding HI supershell can exceed grav- 
ity by orders of magnitude (Fig [5)l , for reasonable assump- 
tions about the gravitational force. It is therefore possible 
that Lya radiation pressure plays an important role in de- 
termining the kinematics of HI supershells around starburst 
galaxies. We have demonstrated that the total Lya radia- 
tion force on a spherical HI supershell can be written as 
Fra.d = MpLa/c, where the 'force-multiplication factor' Mf 
relates to the average trapping time of Lya photons in the 
neutral medium. The factor Mf can greatly exceed unity, 
as illustrated by Fig [6] For comparison, the maximum pos- 
sible radiation force due to continuum radiatiorjj is Lboi/c, 
in which Lboi is the bolometric luminosity of the central 
galaxy. For a typical star fo rming galaxy. La ~ 0.07I/boi (e.g. 
IPartridge fc Peebleslll967^ ■ whereas for a galaxy that con- 
tains population III stars. La ~ 0.24Lboi (Schaerer. 2003). 
Hence, the Lya radiation pressure can dominate the maxi- 
mum possible continuum radiation pressure if Mf >14 (for 
a normal stellar population) , a threshold which is easily ex- 
ceeded at large column densities of relatively slow-moving 
HI shells (see FiggJ. 

The possibility that Lya radiation alone can result 
in a radiation forc e that exceeds Lboi/c is important. 
[Murray et al.l (|2005l ') have shown that the total momentum 
carried by radiation from a star-forming region can exceed 
the total momentum deposited by supernova explosions in 
it, and so galactic outflows may be driven predominantly 
by continuum radiation pressure. We have argued that Lya 
radiation pressure may in some cases be even more impor- 
tant than continuum radiation pressure, and thus provide 
the dominant source of pressure on neutral hydrogen in the 
ISM. 

The important implication of our last result is that Lya 
radiation pressure may drive outflows of HI gas in the ISM. 
Observations of local starburst galaxies have shown that the 
presence of outflowing HI gas may be required to avoid com- 
plete destruction of the Lya radiation by dust and to al- 
low its e scape fro m the host galaxi e s (IKu nth et al. 19981: 
iHaves et a l. 2008,; lOstUn et all l2008l : IXtek et al. ,2008). At 
high redshifts, the Lya emission line of galaxies is often red- 
shifted relative to other nebular recom bination lines (such 
as Ha) and metal a bsorption lines (e.g. iPettini et al.ll200ll : 
IShaplev et al.ll2003l ). Furthermore, the spectral shape of the 



^ Continuum radiation may exert a force on steady-state out- 
flows (with a constant mass ejection rate, M) around late-type 
stars that may significantly exceed MAi) 2> ^bol/c llSalpeten 
ll97'4lMlvezic &: Elitzur 1995J). This is not because of 'trapping' of 
continuum photons, but related to the propagation speed of the 
photons and the wind. We similarly expect the radiative force 
of trapped Lya photons in steady-state outflows to potentially 
exceed MAt) 2> MFLa/c (and as argued in this paper, it is pos- 
sible that MpLa/c > Lbol/c). Note though that these steady- 
state outflows are clearly different from those discussed in § 13.31 
in which a well deflned thin shell HI gas is physically separated 
from the central Lya source. 



Lya emission line is typically asymmetric, wi th emission ex- 
tend ing well into the red wing of the line (e.g. iLegueux et al.l 
|l995). Both of these observations can be explained simul- 
taneously if the observed Lya photons scatter off neu- 
tral hydrogen atoms in an outflowing 'su pershell' of HI 
gas that surrounds the star forming region s ('Lcgucu x et al.l 
1995; Tcnorio-Tag le et al.lll999l : lAhn et aLi r2003: Ahnl l2004l : 
Iverhamme et al.ll2006l , l2008 fl. The possibility that Lya radi- 
ation pressure may be important in determining the proper- 
ties of expanding supershells is exciting, and is discussed in 
more detail in a companion paper (Dijkstra & Loeb 2008). 
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APPENDIX A: HI LEVEL POPULATIONS 

Al The (De)Excitation Rates of the 2p Level 

Our calculations assumed that all neutral hydrogen atoms 
populate their electronic ground state (i.e. the Is level). Be- 
low, we explore the physical conditions under which this 
assumption holds. 

Processes that populate the 2p level include: (i) coUi- 
sional excitation of neutral hydrogen atoms in both the Is 
and 2s states by electrons and protons, (ii) recombination 
into the 2p state following photoionization or coUisional ion- 
ization, (Hi) photoexcitation of 2p level, {iv) photoexcitation 
of np (n > 3) levels, followed by a radiative cascade that 
passes through the 2p state. Processes that de-populate the 
2p level include (v) coUisional de-excitation to by electrons, 
(vi) stimulated Lya emission 2p -^ Is following the absorp- 
tion of a Lya photon, and (vii) spontaneous emission of a 
Lya photon. 

Quantitatively, the rate at which the population in the 
2p level is populated is given by 



dn2p 
dt 



= Cls2p".e".ls + C2s2pnpn2s + 0.68nenHIiarcc,B + (Al) 



2_^ Prinis fn2 — C2plsnen2p — C2p2sWpn2p — P2".2p — n2pA21, 
n=2 

where Ci„'s denote coUisional excitation rate coefficients 
from level I to u (in cm'' s~^), n^ denote number densities 
is species 'x' (i.e rie denotes the electron number density, 
while TiHii denotes the number density of HH ions) , P„ de- 
note photoexcitation rates to the state np (in s~^), and fn2 
denotes the probability that photoexcitation of the np level 
results in a radiative cascade that passes through the 2p 
level. 

In the reminder of this Appendix, we will estimate the 
order-of-magnitude of each term. 

• The Einstein-A coefficient of the Lya transition is 
A21 = 6.25 X 10* s"'^. That is, spontaneous emission of Lya 
depopulates the 2p state at a rate A21 = 6.25 x 10** s~^. 

• CoUisional excitation of neutral hydrogen atoms in Is 
level by electrons populates the 2p level at a rate Cis2p»^e = 

ne[ * -'^^r/2°'' j [^^^^fir] e-^/"^ s-i (e.g. lOsterbrockl 
Il989f ). Here, T denotes the gas temperature in K, rie is the 
electron density in cm"'', gris = 1 and g2p = 3 are the statis- 
tical weights of the Is and 2p levels, flfls, 2p ) = 0.40 — 0.50 
(5000 K < T < 2x10'' K. IOsterbroc^il989f ). and x = 10.2 
eV is the energy difference between the Is and 2p levels. For 
temperatures T < 2x 10* K, we find that Cis2pn-e < 10~^°ne 



• The coUisional excitation rate of neutral hydrogen 
atoms in 2s level by protons (which dominate over collisions 
with electrons by about a factor of ~ 10) populates the 2p 
level at a rate C2. 



2pnp 



2 X 10~^np s"^ (e.g. Oste rbrocM 
[l989|). To assess the term C2s2pfipn-2s requires one to com- 
pute n2s- The fraction of atoms in the 2s state is deter- 
mined by rates similar to those mentioned above, except 
that the 2s-state is metastable and its Einstein coefficient 
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is A2sis ~ 8 s~^. The 2s-state may therefore be overpopu- 
lated relative to the 2 p state by orders of magnitude (see e.g. 
IPennison et al.ll2005l . and references therein). In close prox- 
imity to a luminous source, the 2s level is populated mostly 

via transitions of th e form Is — > 3p -? 2s (jSethi et al.ll2007l : 
iDiikstra et~ai1l2008l '). and n2s ~ 0.l2nisP3/A2si8, where P3 
is the rate at which Ly/3 photons are scattered and the pref- 
actor 0.12 denotes the probability that absorp tion of the Ly/3 
is fol lowed by re-emission of an Ha photon (jDiikstra et al.l 
I2OO8I ). 

• The recombination rate into the 2p level is given by 
0.68 ai.ec.Bnenp = 18 x 1 0-^^(rgaB/10^ Ky^-'^n^rij, cm^ s'^ 
(e.g. lHui fc Gnedinlll997l ). where Up is the proton density in 
cm"^. 

• The rate at which transitions of the form Is -^ np occur 
by absorbing a photon is given by P„ = 47r J ^-^a„{u)du. 
Assuming for simplicity that J(i') does not vary with fre- 
quency, i.e. J(i^) — J, we have 



47rJ/n7re^ 

hPnirieC 



(A2) 



where /„ denotes the oscillator strength of the transition, 
e (nie) the charge (mass) of the electron, and hun denotes 
the energy difference between the Is and np levels. The os- 
cillator strengt h decreases rapidly with i ncreasing n (e.g. 
chapter 10.5 of lRvbicki fc Lightmanlll979l l. and in practice 
we can safely ignore all terms with n > 2. We then need not 
worry about the factors fn2 (which have b een computed by 
iPritchard fc Furlanettoll2006l : lHiratall2006l ). 
The Lya scattering rate is given by 



4:TTr^ Avhiyaniec' 



(A3) 



where we replaced J (in erg s~^ Hz~^ sr~^ cm~^) with 
J = Mf iK^i°A^^ , in which La is the Lya luminosity of the 
central source (in erg s~^), and Au is the frequency range 
over which these Lya photons have been emitted. The fac- 
tor Mf takes into account the fact that resonant scattering 
traps Lya photons in an optically thick medium (see ij |3.3|) . 
Substituting fiducial numbers 



Pa = 2 X 10^ s"^ 



La 



10*2 erg/s 



w10_^Wpc\2/Mf\ 
)\ Au )\r ) \ 



100/ 
(A4) 

The rate at which Ly/3 photons scatter can be related to 
the Lya luminosity, L^, and the equivalent width (EW) of 
the line, if one writes the specific intensity Jivij) near the 
Ly/3 resonance in terms of the Lya luminosity of the central 
source as Jif/s) = -^ ^^ -, g \ i (note that we assumed that 
the specific intensity of the continuum remains constant be- 
tween Va and vp). The rate at which Ly/3 photons scatter 
can then be written as 



2 X 10" 



-'( 



La 



1042 erg s 



.)( 



EW \-^( r \ 
200 A/ Vpc/ 



(A5) 

Equation (|A5|) illustrates that it is very difficult to bring the 
ratio n2s/nis ~ 0.12P3/(A2sis) to unity. 

• CoUisional deexcitation rates relate to th e coUisional 
excitation rates via C;„ — C„;-2^e~-^''°"^ (e.g. lOsterbrockl 
1 19891 ). Combined with the formulas given above, it is 
straightforward to verify that the coUisional de-excitation 



rates are subdominant relative to the rate at which sponta- 
neous Lya emission de-populates the 2p level. 

• Lastly, Eq HA4|) shows that the stimulated emission rate 
isP2 < A21. 



A2 HI Level Populations in this Paper 

The maximum number density of hydrogen nuclei in this 
paper is encountered in § 13.31 for the expanding shell of HI 
gas with TVhi ~ 10^^ cm"^ and a thickness dr — 0.01 kpc, in 
which Timax ~ 30 cm~^. At these column densities, the shell 
self-shields against ionizing radiation, and is likely mostly 
neutral. For simplicity, let us assume that n^ = Up — nni ~ 



30 cm . Under these conditions: 

• the coUisional excitation rate from Is 
-1 



2p is Ci2 < 



• the coUisional excitation rate from 2s — *■ 2p per atom 
in the Is state- is C2s2pn2snp/nis ~ 5 x 10~^n2s/'nis 
s~^— 5 X 10~^P3/A2ais s~^. The maximum Lya luminos- 
ity considered in this paper is ~ 10'*'' erg s~^. Let us 
conservatively assume that EW= 20 A(rest-frame), which 
corresponds roughly to the detection threshold tha t exists 
in narrow-band surveys (e.g. IShimasaku et al.ll2006l ). Using 
EalXSl we find that Ps.max ~ lO"'^ s"^ (for r = 0.1 kpc), 
and therefore that C2s2pn2snp/n\s ~ 10~* s"^. 

• the recombination rate is 5 X 10"^^ (Tgas/lO"* K)"°''s"\ 

• the maximum Lya scattering rate is (substituting r — 
0.1 kpc, Mf = 100, Au = 0.00l!/„ into EqUU P2,max = 0.2 



By comparing these rates to the rate at which spon- 
taneous emission of Lya depopulates the 2p state, A21 — 
6.25 X 10^ s"'^, we find that all excitation rates are ^ 9 
orders of magnitude smaller than the de-excitation rate for 
the wind models discussed in i; l3.3l In equilibrium, hydrogen 
atoms in their electronic ground (Is) state are therefore ^ 9 
orders of magnitude more abundant than hydrogen atoms in 
their first excited (2p) state. Furthermore, as was mentioned 
above the ratio of atoms in the 2s and Is levels is given by 
n2a/nis — 0.12P3,niax/j42sis ~ 10~^. Slucc the densities, 
Lya luminosities, and the Ly/3 scattering rates are lower, 
the fraction of HI atoms in their first excited states are even 
smaller in other sections of the paper. In conclusion, for all 
applications presented in this paper, no accuracy is lost by 
assuming that all hydrogen atoms occupy their electronic 
ground state. 

Finally, in l2.2l we computed solutions for the radial de- 
pendence of U{r) (Fig [T|. The energy density, U{r), was 
quoted to depend linearly on the luminosity of the cen- 
tral source. This is valid unless (i) the Lya scattering rate, 
P2 > ^21 = 6.25 X 10* s"\ or (ii) the Ly/3 scattering 
rate exceeds P3 ^10^ s~^. In either case, our assumption 
that all hydrogen atoms populate their electronic ground 
state breaks down. Condition (i) translates to La ^3 x lO''** 

erg s"^(Ai//0.001i/<,)(r/pc)2(100/AfF) (Eq.ED), while con- 
dition (m) translates to La ^lO"® erg s"^(EW/20 k){r /pcf 
(Eq.lXSl. Substituting r = 0.1 kpc, Mf = 100 (Fig [f] shows 
that resonant scattering enhances the energy density by a 
factor of ^10 relative to La/^Tvr'^c), Av = 10~''i'Q (ther- 
mal broadening alone in T = lO'* K gas results in Av ^ 26 
km s~^, which translates to Av — lO""*!^^), and EW — 20A, 
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condition (i) translates to La >3 x 10^^ erg s^^, while condi- 
tion (ti) translates to La ^W'''^ erg s^^ The more conserva- 
tive condition (ii), La ^10^" erg s~^, translates to N^,4 i^lO^, 
well beyond the regime considered in this paper. 



